Buffered Coordination Modulation as a Means of Controlling Crystal Morphology and Molecular Diffusion in an Anisotropic Metal-Organic Framework.

Significant advances have been made in the synthesis of chemically selective environments within metal-organic frameworks, yet materials development and industrial implementation have been hindered by the inability to predictively control crystallite size and shape. One common strategy to control crystal growth is the inclusion of coordination modulators, which are molecular species designed to compete with the linker for metal coordination during synthesis. However, these modulators can simultaneously alter the pH of the reaction solution, an effect that can also significantly influence crystal morphology. Herein, noncoordinating buffers are used to independently control reaction pH during metal-organic framework synthesis, enabling direct interrogation of the role of the coordinating species on crystal growth. We demonstrate the efficacy of this strategy in the synthesis of low-dispersity single-crystals of the framework Co2(dobdc) (dobdc4-= 2,5-dioxido-1,4-benzenedicarboxylate) in a pH 7-buffered solution using cobalt(II) acetate as the metal source. Density functional theory calculations reveal that acetate competitively binds to Co during crystallization, and by using a series of cobalt(II) salts with carboxylate anions of varying coordination strength, it is possible to control crystal growth along the c-direction. Finally, we use zero length column chromatography to show that crystal morphology has a direct impact on guest diffusional path length for the industrially important hydrocarbon m-xylene. Together, these results provide molecular-level insight into the use of modulators in governing crystallite morphology and a powerful strategy for the control of molecular diffusion rates within metal-organic frameworks.

Preserving Porosity of Mesoporous Metal-Organic Frameworks through the Introduction of Polymer Guests.

High internal surface areas, an asset that is highly sought after in material design, has brought metal-organic frameworks (MOFs) to the forefront of materials research. In fact, a major focus in the field is on creating innovative ways to maximize MOF surface areas. Despite this, large-pore MOFs, particularly those with mesopores, continue to face problems with pore collapse upon activation. Herein, we demonstrate an easy method to inhibit this problem via the introduction of small quantities of polymer. For several mesoporous, isostructural MOFs, known as M2(NDISA) (where M = Ni2+, Co2+, Mg2+, or Zn2+), the accessible surface areas are increased dramatically, from 5 to 50 times, as the polymer effectively pins the MOFs open. Postpolymerization, the high surface areas and crystallinity are now readily maintained after heating the materials to 150 °C under vacuum. These activation conditions, which could not previously be attained due to pore collapse, also provide accessibility to high densities of open metal coordination sites. Molecular simulations are used to provide insight into the origin of instability of the M2(NDISA) series and to propose a potential mechanism for how the polymers immobilize the linkers, improving framework stability. Last, we demonstrate that the resulting MOF-polymer composites, referred to as M2(NDISA)-PDA, offer a perfect platform for the appendage/immobilization of small nanocrystals inside rendering high-performance catalysts. After decorating one of the composites with Pd (average size: 2 nm) nanocrystals, the material shows outstanding catalytic activity for Suzuki-Miyaura cross-coupling reactions.

Ab Initio Flexible Force Field for Metal-Organic Frameworks Using Dummy Model Coordination Bonds.

We present force fields developed from periodic density functional theory (DFT) calculations that can be used in classical molecular simulations to model M-MOF-74 (M = Co, Fe, Mg, Mn, Ni, Zn) and its extended linker analogs. Our force fields are based on cationic dummy models (CDMs). These dummy models simplify the methodology required to tune the parameters and improve the accuracy of the force fields. We used our force fields to compare mechanical properties across the M-MOF-74 series and determine that increasing the size of the linker decreases the framework rigidity. In addition, we applied our force fields to an extended linker analog of Mg-MOF-74 and characterized the free energy of a previously reported deformation pattern in which the one-dimensional hexagonal channels of the framework become irregular. The free energy profiles confirm that the deformation is adsorbate induced and impossible to access solely by a pressure stimulus. On the basis of our results, we conclude that the force fields presented here and others that may be developed using our methodology are transferable across metal-organic framework series that share a metal center topology. Finally, we believe that these force fields have the potential to be adapted for the study of complex problems in MOF chemistry, including defects and crystal growth, that have thus far been beyond the scope of classical molecular simulations.

An experimental and computational study of CO2 adsorption in the sodalite-type M-BTT (M = Cr, Mn, Fe, Cu) metal-organic frameworks featuring open metal sites.

We present a comprehensive investigation of the CO2 adsorption properties of an isostructural series of metal-organic frameworks, M-BTT (M = Cr, Mn, Fe, Cu; BTT3- = 1,3,5-benzenetristetrazolate), which exhibit a high density of open metal sites capable of polarizing and binding guest molecules. Coupling gas adsorption measurements with in situ neutron and X-ray diffraction experiments provides molecular-level insight into the adsorption process and enables rationalization of the observed adsorption isotherms. In particular, structural data confirms that the high initial isosteric heats of CO2 adsorption for the series are directly correlated with the presence of open metal sites and further reveals the positions and orientations of as many as three additional adsorption sites. Density functional theory calculations that include van der Waals dispersion corrections quantitatively support the observed structural features associated with the primary and secondary CO2 binding sites, including CO2 positions and orientations, as well as the experimentally determined isosteric heats of CO2 adsorption.

Unexpected Diffusion Anisotropy of Carbon Dioxide in the Metal-Organic Framework Zn2(dobpdc).

Metal-organic frameworks are promising materials for energy-efficient gas separations, but little is known about the diffusion of adsorbates in materials featuring one-dimensional porosity at the nanoscale. An understanding of the interplay between framework structure and gas diffusion is crucial for the practical application of these materials as adsorbents or in mixed-matrix membranes, since the rate of gas diffusion within the adsorbent pores impacts the required size (and therefore cost) of the adsorbent column or membrane. Here, we investigate the diffusion of CO2 within the pores of Zn2(dobpdc) (dobpdc4- = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) using pulsed field gradient (PFG) nuclear magnetic resonance (NMR) spectroscopy and molecular dynamics (MD) simulations. The residual chemical shift anisotropy for pore-confined CO2 allows PFG NMR measurements of self-diffusion in different crystallographic directions, and our analysis of the entire NMR line shape as a function of the applied field gradient provides a precise determination of the self-diffusion coefficients. In addition to observing CO2 diffusion through the channels parallel to the crystallographic c axis (self-diffusion coefficient D∥ = (5.8 ± 0.1) × 10-9 m2 s-1 at a pressure of 625 mbar CO2), we unexpectedly find that CO2 is also able to diffuse between the hexagonal channels in the crystallographic ab plane (D⊥ = (1.9 ± 0.2) × 10-10 m2 s-1), despite the walls of these channels appearing impermeable by single-crystal X-ray crystallography and flexible lattice MD simulations. Observation of such unexpected diffusion in the ab plane suggests the presence of defects that enable effective multidimensional CO2 transport in a metal-organic framework with nominally one-dimensional porosity.

The Influence of Intrinsic Framework Flexibility on Adsorption in Nanoporous Materials.

For applications of metal-organic frameworks (MOFs) such as gas storage and separation, flexibility is often seen as a parameter that can tune material performance. In this work we aim to determine the optimal flexibility for the shape selective separation of similarly sized molecules (e.g., Xe/Kr mixtures). To obtain systematic insight into how the flexibility impacts this type of separation, we develop a simple analytical model that predicts a material's Henry regime adsorption and selectivity as a function of flexibility. We elucidate the complex dependence of selectivity on a framework's intrinsic flexibility whereby performance is either improved or reduced with increasing flexibility, depending on the material's pore size characteristics. However, the selectivity of a material with the pore size and chemistry that already maximizes selectivity in the rigid approximation is continuously diminished with increasing flexibility, demonstrating that the globally optimal separation exists within an entirely rigid pore. Molecular simulations show that our simple model predicts performance trends that are observed when screening the adsorption behavior of flexible MOFs. These flexible simulations provide better agreement with experimental adsorption data in a high-performance material that is not captured when modeling this framework as rigid, an approximation typically made in high-throughput screening studies. We conclude that, for shape selective adsorption applications, the globally optimal material will have the optimal pore size/chemistry and minimal intrinsic flexibility even though other nonoptimal materials' selectivity can actually be improved by flexibility. Equally important, we find that flexible simulations can be critical for correctly modeling adsorption in these types of systems.

Adsorbate-induced lattice deformation in IRMOF-74 series.

IRMOF-74 analogues are among the most widely studied metal-organic frameworks (MOFs) for adsorption applications because of their one-dimensional channels and high metal density. Most studies involving the IRMOF-74 series assume that the crystal lattice is rigid. This assumption guides the interpretation of experimental data, as changes in the crystal symmetry have so far been ignored as a possibility in the literature. Here, we report a deformation pattern, induced by the adsorption of argon, for IRMOF-74-V. This work has two main implications. First, we use molecular simulations to demonstrate that the IRMOF-74 series undergoes a deformation that is similar to the mechanism behind breathing MOFs, but is unique because the deformation pattern extends beyond a single unit cell of the original structure. Second, we provide an alternative interpretation of experimental small-angle X-ray scattering profiles of these systems, which changes how we view the fundamentals of adsorption in this MOF series.

Acetylation within the First 17 Residues of Huntingtin Exon 1 Alters Aggregation and Lipid Binding.

Huntington's disease (HD) is a genetic neurodegenerative disorder caused by an expanded polyglutamine (polyQ) domain near the N-terminus of the huntingtin (htt) protein. Expanded polyQ leads to htt aggregation. The first 17 amino acids (Nt(17)) in htt comprise a lipid-binding domain that undergoes a number of posttranslational modifications that can modulate htt toxicity and subcellular localization. As there are three lysines within Nt(17), we evaluated the impact of lysine acetylation on htt aggregation in solution and on model lipid bilayers. Acetylation of htt-exon1(51Q) and synthetic truncated htt-exon 1 mimicking peptides (Nt(17)-Q35-P10-KK) was achieved using a selective covalent label, sulfo-N-hydroxysuccinimide (NHSA). With this treatment, all three lysine residues (K6, K9, and K15) in Nt(17) were significantly acetylated. N-terminal htt acetylation retarded fibril formation in solution and promoted the formation of larger globular aggregates. Acetylated htt also bound lipid membranes and disrupted the lipid bilayer morphology less aggressively compared with the wild-type. Computational studies provided mechanistic insights into how acetylation alters the interaction of Nt(17) with lipid membranes. Our results highlight that N-terminal acetylation influences the aggregation of htt and its interaction with lipid bilayers.

Communication: Relationship between solute localization and diffusion in a dynamically constrained polymer system.

We investigate the link between dynamic localization, characterized by the Debye-Waller factor, 〈u(2)〉, and solute self-diffusivity, D, in a polymer system using atomistic molecular dynamics simulations and vapor sorption experiments. We find a linear relationship between lnD and 1/〈u(2)〉 over more than four decades of D, encompassing most of the glass formation regime. The observed linearity is consistent with the Langevin dynamics in a periodically varying potential field and may offer a means to rapidly assess diffusion based on the characterization of dynamic localization.

The effects of flanking sequences in the interaction of polyglutamine peptides with a membrane bilayer.

Huntington's disease (HD) is caused by the presence of an extended polyglutamine (polyQ) region at the N-terminus of the huntingtin (htt) protein. The presence of flanking sequences adjacent to the polyQ region has been reported to modulate the effects of potentially toxic protein-membrane interactions. In this study, we consider four peptide systems with various combinations of flanking sequences (KKQ35KK, KKQ35P11KK, N17Q35KK, N17Q35P11KK) and use atomistic molecular dynamics simulations to study the interactions with a DOPC lipid bilayer. We observe significant membrane thinning, disorderliness of lipid molecules, and compensation effects between the top and the bottom leaflets of the bilayer depending on the presence of particular flanking sequences. Overall, we find that the presence of the N-17 flanking sequence is crucial for membrane interactions. Polyproline decreases the interaction with the membrane in the absence of N-17, but enhances it when present along N-17.